Pattern collapse as a mechanism for the formation of solitary structures 
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We report a new mechanism for the formation of localized states, which takes place without front 
propagation. Correspondingly, localized structures appear as solitary states, displaying a behavior 
of single independent cells. The phenomenon is observed in the liquid crystal light-valve experiment 
and is described by a one-dimensional normal form model. We show that such solitary structures 
exist when a pattern solution collapses and its ghost remains to influence the phase portrait. 



The concept of soliton as a single coherent structure 
originates from Hamiltonian systems [l[ . The generaliza- 
tion of this concept to dissipative and out of equilibrium 
systems has led to several studies in the last decades and 
to the definition of a class of localized structures that are 
intended as patterns appearing in a restricted region of 
space. Localized structures have been observed in dif- 
ferent fields, such as in magnetic materials, chemical re- 
actions, granular media, plasmas, nonlinear optics and 
surface waves 0. Theoretically, they are understood 
as macroscopic particle-like objects arising in the pres- 
ence of metastability between different states. Localized 
structures, thus, realize the spatial connection between 
two possible states and appear in the pinning range of 
the front solution Q. In one-dimensional systems and 
from the geometrical point of view, they are described 
by homoclinic orbits passing close to a pattern state and 
converging to an homogeneous state d,Q. More recently, 
they have been explained in terms of front interactions 
@ and their existence has been generalized to the case 
when the homoclinic orbits connect two different pattern 
states, thus leading to localized peaks 0]. 

In all these theoretical frameworks, localized structures 
are predicted to exist in parameter regions where the sys- 
tem exhibits coexistence between two spatially extended 
states and are expected to appear (disappear) by a se- 
quence of saddle- node bifurcations, each bifurcation lead- 
ing to a larger (smaller) number of pattern cells and oc- 
curring in the pinning range of the Pomeau front [3[ . Nev- 
ertheless, there is a wealth of experimental observations, 
for example in optics [8|, [9( or in vibrated granular me- 
dia [10( , of solitary states that display more a character 
of single objects than that of a class of differently sized 
patterns, that is, only the single-cell localized structures 
appear. Recent numerical simulations have also shown 
that cavity solitons occurring in semiconductor micro- 
cavities have their region of existence outside the range 
of bistability between the homogeneous and the pattern 
states , which is consistent with the experimental ob- 
servations ■ The common property of all these obser- 
vations is that localized structures seem to exist without 
the need of metastability, which claims for searching new 
mechanisms for their formation. 

The aim of this letter is to show a new type of transi- 



tion from extended patterns to localized structures, tak- 
ing place without front propagation and leading to soli- 
tary structures that appear randomly in space, behave 
as independent particles and interact through the oscil- 
lations on their tails. We explain this transition through 
the collapse of the pattern solution via a saddle-node bi- 
furcation, which implies the destruction of the homoclinic 
snaking sequence associated to the whole family of dif- 
ferently sized localized patterns. After the collapse, the 
ghost of the pattern remains to influence the phase por- 
trait and solitary (one-cell) structures exist as the result 
of the unique (single- loop) surviving homoclinic orbit. 
We present experimental evidence of this transition for a 
liquid crystal light-valve (LCLV) with optical feedback. 
The results are confirmed by numerical simulations of 
the full model. Then, we introduce a ID normal form 
equation, which allows the use of general geometrical ar- 
guments to describe the pattern collapse. 
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FIG. 1: Phase diagram Ii„ vs Vo- Continuos lines correspond 
to the numerically calculated boundaries from crystal to soli- 
tary structures and to homogeneous state. The dashed line 
delimits the region of bistability between the homogeneous 
and the crystal state. Dots are experimental points. 

The experiment consists of a LCLV inserted in an op- 
tical feedback loop [l3[ . When a voltage Vb is applied to 
the valve, the voltage that effectively drops across the 
liquid crystals is Vlc = rVb + al w , where I w is the 
light intensity impinging on the photoconductor side of 
the LCLV and T, a are phenomenological parameters 
summarizing, in the linear approximation, the response 
of the photoconductor. The input beam, when passing 
through the liquid crystal layer, undergoes a phase shift 
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ip — /3cos 2 #, with 9 the average tilt angle of the liquid 
crystal molecules and (3 = 2irAnd/A, where A is the op- 
tical wavelength and An = n e — uq is the liquid crystal 
birefringence, n e and n being, respectively, the extraor- 
dinary (parallel to the liquid crystal director ft) and ordi- 
nary (perpendicular to ft) refractive index. The tilt angle 
9 obeys a Debye relaxation equation [1J| , 



rdtO = l 2 V\9 -9 + f(0) 



(1) 



where I = 20 /im is the diffusion length summarizing 
the elastic coupling in the liquid crystal and the charge 
diffusion in the photoconductor, r = iO ms is the local 
relaxation time, Vft = 1.05 V is the Freedericksz transi- 
tion voltage and f(9) = n/2 (l — \/Vft /Vlc^J accounts 

for the response of the LCLV [l5|. The optical feedback 
is obtained by sending back onto the photoconductor the 
light that has passed through the liquid crystals and has 
been reflected by the LCLV. After free propagation and 
polarization interference, the light intensity arriving at 
the photoconductor is 

Ii 



III) — 



(2) 



where Jj n is the input intensity and L is the free propa- 
gation length in the optical feedback loop. 
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FIG. 2: Experimental PDF (left) and snapshots (right) show- 
ing the transition from a) an extended pattern to c) a frozen 
configuration of solitary structures, through b) their interac- 
tion and annihilation. 

In the experiment, we have fixed the free propagation 
length to L — 8 cm and we change the applied voltage Vq 



and the input light intensity as control parameters. 
By varying either Vq or Ii n we observe different regimes 
of patterns and localized structures. The transition from 
an extended crystal-like pattern to fioccur in a wide range 
of the Vq, Iin parameters. FigfT] shows the region of ex- 
istence of Bin the Vq — U n plane. Dots are experimental 
points obtained by detecting the pattern change when 
varying the voltage Vo for each value of the pump inten- 
sity Ii n ■ Lines are numerically calculated from the model 
equations Eqs.|T|) and $Z$i, in the following way. We have 
fixed the parameters in the region of existence of soli- 
tary localized structures (SLS), then a single structure 
is switched on through a light pulse. The voltage is in- 
creased until the transition toward hexagons occurs or is 
decreased until the homogeneous steady state (HSS) is 
reached. The same procedure is repeated starting from 
the extended pattern and by decreasing the voltage un- 
til the homogeneous solution is approched. The contin- 
uos lines marks the boundaries of the SLS region. The 
dashed line delimits the region of bistability between the 
homogeneous and the pattern state. This region is very 
thin and is not observed in the experiment due to liquid 
crystal inhomogeneities. Note also that for input inten- 
sity Iin > 0.8 mW/cm 2 a transition from hexagons to 
space-time chaos occurs and for voltage Vq < 13.05 V we 
observe the formation of triangular localized structures 
[l6| . In this case the transition to extended pattern shows 
a more complex scenario, that will not be reported here. 

Experimentally, the boundary between HSS and SLS 
is found by applying a local perturbation at the HSS 
state and by checking if one or more solitary structures 
remain after removing the perturbation, while the bound- 
ary between the crystal and the SLS is determined by 
starting with Vq in the crystal region and then decreas- 
ing it to the SLS region. We show in Fig[2] three typi- 
cal snapshots recorded for = 0.38 mW/cm 2 and by 
changing the rms value of Vo from 13.30 to 13.22 V . 
We first observe an hexagonal pattern (FigJ^h)- When 
decreasing the voltage, a transition takes place to a fi- 
nal distribution of solitary structures appearing in ran- 
dom space positions (Fig[2t). During the transition soli- 
tary single-cell structures interact one with each other 
through the low amplitude oscillations around their pro- 
file (Figl^b). This state displays a gas- like behavior, in 
the sense that it is mainly characterized by continuous 
collisions between the particles. However, at difference 
with a real gas, the collisions may lead also to particle 
annihilation, so that during the time the total number 
of particles is not conserved and their mean distance in- 
creases. After a transient, which lasts for a few seconds, 
a final frozen configuration of particles is reached, where 
firemain fixed in their position. Starting from different 
initial conditions different final frozen configurations are 
observed, so that we can identify the state of FigJ5b as one 
with a large configurational entropy, in the sense defined 
in Ref . [l7j ■ Note that fialways remain individuals and do 
not form a pattern, as can be remarked by the empty 
space between them. We can also note in FigJ^b (at 
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difference with Fig[2^) the presence of oscillatory rings 
around each individual cell, these rings being due to the 
diffraction of a single spot over the uniform background 
and, thus, representing a strong signature of the absence 
of front propagation during the transition. 
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FIG. 3: Experimental profiles recorded during a two-particle 
collision: after the collision one particle has been annihilated. 



In order to characterize the different states, we have 
recorded long-time movies and from each movie we have 
extracted the probability distribution function (PDF) of 
the light intensity fluctuations SI = I— < I >, where I 
is the intensity of each pixel and < I > is the average 
intensity distribution, which is calculated by averaging 
pixel by pixel over the entire stack of images. The re- 
sults are shown in the left part of FigO The crystal and 
the frozen configuration of Bare almost stationary states 
and the PDFs are given by fluctuations of the particle 
around their equilibrium positions. In the crystal case 
the PDF is larger because of slow sliding and gliding of 
differently oriented domains, whereas for the frozen gas 
only the fluctuations due to inhomogeneities and noise 
naturally present in the system are influencing the sta- 
tionary state. The transition from hexagons to solitary 
structures is characterized by an abrupt change of the 
PDF that becomes strongly asymmetric acquiring a high 
tail at large intensity fluctuations (Figf2b). These large 
fluctuations correspond to events of particle interactions, 
when high intensity pulses are produced during the col- 
lisions of two particles. In Fig [3] is displayed a set of 
experimental profiles recorded during an event of two- 
particle collision. It can be seen that a high intensity 
pulse is produced when the two particles collide and that 
after the collision one particle has been annihilated. 

Fig[4] shows numerical snapshots of liquid crystal tilt 
angle 0(x,y,t) during the transition 1 — ► 2 — * 1, as 
marked in Fig[TJ In FigQJi a single cell structure is 
switched on by a triggering intensity pulse (state 1 in 
FigHJ. Then, the voltage is increased to reach the crys- 
talline structure of FigfJ^ (state 2). Figs0b,c,d show the 
transient towards the crystal state. During this tran- 
sition we observe that new cells nucleate spontaneously 



around the starting one. In the experiment this effect 
is not observed because of the presence of spatial inho- 
mogeneities in the LCLV, thus the structures are spon- 
taneously created in random space positions. When the 
voltage is decreased back to its initial value, state 1, an 
ensemble of fiis obtained, as displayed in FigUf. The 
transition from crystal to fidisplays the same scenario 
as the one observed in the experiment, with a transient 
characterized by particle interaction. Correspondingly, 
the PDF of the intensity obtained from the numerical 
simulations are in good agreement with the experimen- 
tal PDFs. 




FIG. 4: Numerical snapshots showing the evolution of the 
liquid crystal tilt angle 6(x,y,t) during the transition 1 — ► 
2 — > 1 indicated in Fig[T] for I in = 0.45 mW/cm 2 . a) Vb = 
13.30 V (state 1); b),c),d),e) time-sequence to reach the state 
2, Vb = 13.22 V; f) state 1, same values of Vb as in a). 



To describe a possible mechanism for the formation of 
solitary structures we consider a simple one-dimensional 
model based on a Ginzburg-Landau equation with spa- 
tial forcing Q . The required ingredients of the model are 
the coexistence between an uniform and a pattern state 
and a saddle-node bifurcation leading to the collapse of 
the pattern solution under the variation of a control pa- 
rameter. We take the following normal form 

d t A = eA- v\A\ 2 A + a\A\*A + \A\ 6 A + d xx A + r)A 2 e iqx , 

(3) 

which is a generic one and can be shown to apply to 
the case of the LCLV for the range of parameters corre- 
sponding to the region of observation of solitary struc- 
tures. A(x, t) is a complex amplitude, e is the bifurca- 
tion parameter, (v, a) control the type of the bifurcation, 
r\ is the intensity of the spatial forcing and q is the wave 
number of the spatially periodic forcing. For negative 
e, v, a and small 77, the system has a stable uniform 
state, \A\ = 0, plus one stable and two unstable pattern 
states. Numerically, we have obtained solitary structures 
when the stable and unstable pattern collapse through a 
saddle-node bifurcation. In correspondence, all patterns 
collapse to the smallest (single-cell) one and only soli- 
tary structures are observed, which continue to exist in a 
large region outside the region of bistability between the 
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homogeneous and the pattern state. 

a) 




W°(H) " 



,'W U (P) 



UP 




FIG. 5: Evolution of the homoclinic structure in the phase 
portrait of the Poincare section: a) the homoclinic snaking 
occurring in the pinning range and far from the saddle-node 
bifurcation (77 = 0.3,e = —2.6), b) close (rj = 0.3,e = —2.5) 
and c) after the saddle- node bifurcation [r] = 0.3,e = —2.4). 
(The other parameters are v — 0.7, a = —4.9, q = 12). 

To illustrate the phenomenon of the pattern collapse, 
we show in Figj5] how it is modified the typical sce- 
nario of homoclinic snaking occurring in the Poincare 
section of the spatial reversal stationary dynamical sys- 
tem [l8| . FigJSJi represents the pinning range. IT is the 
symmetry plane related to the spatial reflection symme- 
try (x — > —x), that is, the trajectories are symmetrical 
with respect to this plane. For the model Eq.([3]) the re- 
flection symmetry has the form x — > — x and A — > A, 
the symmetry plane is defined by IT = {$l(d x A) — 
3(A) = 0}. The points H, P, and UP represent, 
respectively, the stable homogeneous solution, the sta- 
ble pattern solution and the unstable pattern solution. 



The continuous and dashed dark {W 8 (H) ,W U (H)} and 
grey {W S (P), W U (P)} curves are the stable and unsta- 
ble manifold of H and P, respectively. The homoclinic 
curves, represented by the points {hi, h%, • • • } correspond 
to the differently sized localized patterns with one, two, 
three cells and so forth, in agreement with the bifurca- 
tion diagram of localized patterns proposed in Q . When, 
by changing a control parameter, UP approaches P, the 
influence of the heteroclinic tangle is decreased, as illus- 
trated in Fig[5b. After the saddle node bifurcation the 
pattern solution does not exists any more and all the ho- 
moclinic solutions collapse in the same one, as depicted in 
Fig. [5}:. The saddle-node bifurcation continues to influ- 
ence the manifold evolution, because the pattern leaves 
its ghost in the phase space, which continues to attracts 
the phase space trajectories [l9j]. As a consequence, the 
remaining homoclinic orbit represents a solitary (single- 
cell) localized structure, which, by topological reasons, 
continues to survive until it collides with another homo- 
clinic solution. Note that a similar mechanism has been 
recently proposed to explain the transition from oscillat- 
ing localized structures to an excitability regime, where 
the collapse takes place through the merging of a limit 
cycle with a saddle point [20j . 

In conclusion, we have shown that pattern collapse is a 
new mechanism for the formation of solitary structures, 
the minimal requirement being the presence of a stable 
solution with a pattern ghost that allows for the existence 
of a single-loop homoclinic orbit. 

The authors thanks the support of ECOS-CONICYT 
collaboration program. M.G. C. acknowledges the finan- 
cial support from the ring program ACT15 of Programa 
Bicentenario. U. Bortolozzo acknowledges a fellowship of 
the Ville de Paris. 



[1] A.C. Newell, Solitons in mathematics and physics, (So- 
ciety for Industrial and Applied Mathematics, Philadel- 
phia, 1985). 

[2] See, e.g. ,M. C. Cross and P. C. Hohenberg, Rev. Mod. 

Phys. 65, 851 (1993) and references therein. 
[3] Y. Pomeau, Physica D 23, 3 (1986). 
[4] P. Coullet, C. Riera, and C. Tresser, Phys. Rev. Lett. 84, 

3069 (2000). 

[5] W. van Saarloos, and P.C. Hohenberg, Phys. Rev. Lett. 

64, 749 (1990). 
[6] M.G. Clerc, and C. Falcon, Physica A 356, 48 (2005). 
[7] U. Bortolozzo, M.G. Clerc, C. Falcon, S. Residori, and 

R. Rojas, Phys. Rev. Lett. 96, 214501 (2006). 
[8] B. Schapers, M. Feldmann, T. Ackemann, and W. Lange, 

Phys. Rev. Lett. 85, 748 (2000). 
[9] P.L. Ramazza, E. Benkler, U. Bortolozzo, S. Boccaletti, 

S. Ducci, and F.T. Arecchi, Phys. Rev. E 65, 066204 

(2002). 

[10] P. B. Umbanhowar, F. Melo, and H. L. Swinney, Nature 
(London) 382, 793 (1996). 



[11] G. Tissoni et al., J. Op. Soc. Am.B. 16, 2095 (1999). 
[12] X. Hachair et al. Phys. Rev. A 69, 043817 (2004). 
[13] S. Residori, Phys. Rep. 416, 201 (2005). 
[14] M. G. Clerc, A. Petrossiam and S. Residori Phys. Rev. 

E 71, 015205 (2005). 
[15] S. Residori, U. Bortolozzo , Mol. Cryst. Liq. Cryst. 

454,207 (2006); U. Bortolozzo, S. Residori, Phys. Rev. 

Lett. 96, 037801 (2006). 
[16] U. Bortolozzo, L. Pastur, P.L. Ramazza, M. Tlidi, G. 

Kozyreff Phys. Rev. Lett 93, 253901 (2004). 
[17] P. Coullet, C. Toniolo, and C. Tresser, Chaos 14, 839 

(2004). 

[18] P.D. Woods, and A.R. Champneys, Physica D 129, 147 
(1999). 

[19] S. H. Strogatz, Nonlinear Dynamics and Chaos: With Ap- 
plications to Physics, Biology, Chemistry and Engineer- 
ing, (Addison- Wesley, Reading, Massachussets,1994). 

[20] D. Gomila, M.A. Matias, and P. Colet, Phys. Rev. Lett. 
94, 063905 (2005). 



